Conformal invariant saturation 



H. Navelet and R. Peschanski* 

We show that, in onium-onium scattering at (very) high energy, a transition to saturation happens 
due to quantum fluctuations of QCD dipoles. This transition starts when the order c? correction of 
the dipole loop is compensated by its faster energy evolution, leading to a negative interference with 
the tree level amplitude. After a derivation of the the one-loop dipole contribution using conformal 
invariance of the elastic 4-gluon amplitude in high energy QCD, we obtain an exact expression of 
the saturation line in the plane (Y,L) where Y is the total rapidity and L, the logarithm of the 
onium scale ratio. It shows universal features implying the Balitskyi - Fadin - Kuraev - Lipatov 
(BFKL) evolution kernel and the square of the QCD triple Pomeron vertex. For large L, only the 
higher BFKL Eigenvalue contributes, leading to a saturation depending on leading log perturbative 
. . . QCD characteristics. For initial onium scales of same order, however, it involves an unlimited 

summation over all conformal BFKL Eigenstates. In all cases, conformal invariance is preserved for 
■ the saturation mechanism based on dipole loops. 
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I. SATURATION FROM QCD DIPOLES 



. A regime of QCD at small coupling constant where the density of partons begins to saturate has since long been the 
' subject of many studies In one class of models, one expects the effect due to the high energy multiplication of 

partons due to the QCD dynamics. Another way which was proposed to investigate high density partonic effects is 
^ the consideration of collisions at high energy on heavy nuclei. Indeed, the number of partons is supposed to be large 
in the background of the collision providing boundary conditions favorable for the saturation mechanism to happen. 
OO In the present paper, we will focus on the first scheme, i.e. the saturation mechanism due to the energy evolution 
of the parton density in a purely perturbative QCD framework. More precisely, we will consider even "perturbative 
QCD" boundary conditions, in order to select the features of saturation which could be fully calculated from the QCD 
' lagrangian. 

On the experimental ground, saturation has not been clearly identified in the evolution of structure functions in the 
HERA range, even if some evidence for partial saturation has been found by analyzing vector meson diffractive pro- 
, duction in impact parameter space [|l9| . The problem of identifying saturation effects in actual or future experiments 
^Y^' is very stimulating, since the increase of the energy range is at reach at Tevatron and later at LHC. Even, some lessons 
O il will be available from experiments with heavy nuclei at RICH. The goal of our purely perturbative calculations, even 
1^ . if they appear to be far from the experimental realizations, is to give an example of saturation which could be fully 
computed, without referring to a non-perturbative ansatz. In practice, later on, it will be possible, using factorization 
properties, to combine our calculations with some reasonable non-perturbative wave functions to get predictions for 
realistic processes, like hard hadron-hadron intcrations at Tevatron and LHC. 

Our idea is to start from the most elementary hard processes, namely onium-onium scattering at (very) high 
' energies. The reaction can be described j|J|] by the wave-functions of the onia in terms of qq colour singlet dipoles 
at the time of the (hard) interaction, which then interact through double gluon exchanges. This description, which 
corresponds to the tree level dipole interactions, is equivalent to the usual Balitsky, Fadin, Kuraev, Lipatov (BFKL) 
formalism pO| for the same process, as can be proven [ pT[ | using the conformal invariant symmetry of the BFKL 
evolution equation. When energy increases, the wave functions has an ever increasing density of colour dipoles and 
the tree level dipole interaction is expected to be modified by unitarity corrections which are already present in the 
per turbative regime. The first order correction is the one-loop dipole contribution which was first considered in Ref. 

as a formal calculation of next-to-next order contributions in the coupling constant. In the present paper we 
show that, despite its higher order, this one-loop unitarity correction gives rise to the onset of saturation in onium- 
onium scattering and thus allows one to give a complete perturbative QCD description of the saturation critical line. 
This saturation mechanism is to be distinguished from the one which would come from the interaction of a onium 
(e.g. a collection of independent dipoles) with a large nucleus (see, e.g. |^,|lT|- |l3| , p^ , p7Hl^ ) , which, in some cases, 
can be associated with a non-linear branching process without dipole recombination. In the case we consider, the 
non-linearities due to the high density of partons correspond to the recombination of QCD dipoles. 
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We write the onium-omum elastic cross-section 



where Q is the transfer momentum. The amphtude reads 



Fq^~i <fpd^p' / dzdz' <^{p,zW,z') fQ{p,p'\Y) , (2) 







where the onium is defined as a heavy qq state with q (resp. q) having the fraction z (resp. 1 — 2;) of the longitudinal 
momentum in an infinite momentum frame. It is defined by its wave-function <&(/9, z) over a basis of dipole states 
of transverse size p = pq — pi, where pq, pi are the positions of the dipole's constituants in transverse coordinate 
space. fQ{p, p'\Y) is the dipole-dipole elastic amplitude for a full rapidity interval Y and a momentum transfer Q, 
the calculation of which at one dipole-loop order and saturation properties are the subject of the present paper. 

In fact the physical origin of the saturation mechanism that we want to investigate is the following: it is known 
that the wave function of an Onium in the infinite momentum frame develop a large number of dipole states as a 
function of the rapidity (which plays the role of time in a cascading mechanism). However, the larger is the dipole 
number, the smaller is their size and thus they do not contribute to saturation. This can be pro ven rigorously by 
the fact that the dipole-dipole cross-section add together to give back the BFKL cross-section It is interesting 
to note the use of conformal invariance to get this exact bootstrap property between an "s-channel" and "t-channel" 
description of the same amplitude. 

However, now if one computes the dipole-dipole interaction at one-loop level, the system is allowed to take into 
account quantum fluctuations where virtual dipole pairs may be created during the time of interaction. Moreover, 
all sizes of dipoles can contribute in the loop, which goes beyond the limitations of the tree-level mechanism. This 
is expected from perturbative unitarity arguments ^ which expect the damping of the energy dependence of the 
amplitude by saturation effects due to multi-loop dipole contributions. This is precisely the goal of our work to put 
this expectation on a quantitative level. 

Through unitarity arguments, the QCD saturation phase corresponding to the mechanism we consider requires an 
estimate of all dipole loop contributions, at least in the high energy approximation. In the present paper we will 
concentrate on the first step, i.e. the transition to saturation which is expected to be due to the interfernece between 
the tree and 1-loop level. We expect that the techniques which we were led to consider in this work, relying notably 
on the conformal invariance properties, will be useful to investigate this infringing phase of matter. 

Our main results are the following: 

i) Up to one-loop QCD dipole level, the dipole-dipole amplitude reads 



fQ{p,p'\Y) ^ + /« = ^ \pp'\ I dh Eq{p')E^{p) d{h) e-C')^ 
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(3) 



In equation (^), the symbolic notation J dh = X]^-oo / '^'^ corresponds to the integration over the quantum 
numbers associated to the continuous unitary irreducible representations of SL{2,C), namely 

1 — n ~ \ -\- n ^ ... 

h = ii' ^ — ; h^ii' -\ — ; n €M , ly eTZ . (4) 

The SL{2,C) Eigenvalue of the BFKL kernel is 

= ^ x{h) ^ ^ {2^{l) -^{h)-^{l-h)}, (5) 

where ^ = (logF)', and Eq{p) are the SL{2,C) Eigcnfunctions gzv{h, ha, hi,) is the known QCD triple Pomeron 
coupling function ||2j,|2^, whose definition is given in Appendix 1, together with Eq{p) and the other ingredients 
(such as the functions a{h) ,d{h)) of the S'L(2, C)-invariant formalism. 

The symmetry of the tree level amplitude remains valid at one-loop level, leading to a conformal invariant saturation 
scheme. 
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ii) The evaluation of formula (g) when the condition L = log 



^ 1 is fulfilled can be done in the forward 



p' 

direction (Q — 0) in terms of the leading Eigenfunction and Eigenvalue of the BFKL kernel, namely with n = 0. By 
contrast, for similar dipole sizes {L — 0(1)) infinitely many Eigenfunctions and Eigenvalues contribute to the loop 
expression, leading to a quite different structure. 

The QCD critical saturation line is derived in the L 3> 1 regime and is given by: 

Y(^^n_^y^^^^^^^^^ ^fc + ^logy-logL, (6) 

where Y = w(l/2)F and ^ C{log(l/a^)} is a large constant (formula ( p^ ) in the text). The transition to 
saturation occurs when the first member of (^ exceeds the second one. Note that using the variable Y means that 
the rapidity (i.e. energy) dependence of the amplitude is quantified by reference to the BFKL intercept 1 + w(l/2). 
Hence the overall F, L dependence is determined by "universal" numbers, independent of the coupling constant. 

In formula (^), there appears two critical values of the effective anomalous dimensions ft.c for the tree level and h* 
for the one-loop contribution. They are both defined as solution of implicit equations 

^ = ^ ; xiK) = 2x(l/2) . (7) 

he is a moving saddle point in the tree-level integrand of (||) while h* ^ .183 is the fixed location of a double pole 
singularity in the one- loop integrand of (^). 

The plan of the paper is the following. The next section II is devoted to the derivation of formula (||) and Section 
III to that of formula (^). A disussion of the resulting saturation line and an outlook are proposed in section IV. The 
S'L(2, C)-invariant formalism which we use is presented in Appendix 1 and some technical parts of the derivation of 
the one-loop amplitude are present in Appendix 2 and 3. 



II. THE ONE-LOOP DIPOLE AMPLITUDE 

p(i) 



At one-loop level, the amplitude /g has a simple and physically appealing formulation [Q in the transverse 
coordinate space, see Fig.l. We write 

fS\p,p'\Y) ^ ^ Jd\ e^Q r(i) {p,p,-p',p',\Y) , (8) 

with impact parameter r. The dipole-dipole amplitude in coordinate space reads 

'' d^Paod^Pa.d^Pbod^Pb, d^Pa'/^Pa'^d^Pb'^d^Pb[ 



T<^'HpoPi;PoP'i\Y = y+y')^- f'dy r dy' f ■ 

Jo Jo J 



\Pa Pb\^ \pa' Pb'l'^ 



X n2{poPl;PaoPa^,PboPb^\y-y,y)n2 {pop[; p'^^p'^^, p'f,^p'i,^\y' -y' ,y') T{paoPai,Pa'„Pa[) T {pb„ Pb, , Pb'gPb[) (9) 

where po, Pi are the transverse coordinates of one of the initially colliding dipoles (resp. p'q ,p\ for the second one), 
while PaoPai and pboPbi, are the two interacting dipoles emerging from the dipole popi after evolution in rapidity 
(resp. P{i'], for the second one). 

The overall minus sign comes from the well-known AGK rules relating the amplitudes to the discontinuity here 
calculated in the c.o.m. frame of the dipole-dipole reaction. The tree level interaction amplitudes T are given in 
Appendix 1. In formula (||), n2 {poPi', PaoPai, PboPbi\y — y,y) stands for the probability of producing two dipoles after 
a mixed rapidity evolution [pp^] , namely with a rapidity y — y with one-Pomeron type of evolution and a rapidity y 
with two-Pomeron type of evolution, see Fig. 2. Further on, one has to integrate over y,y' . The solution for 712 is 
given in Appendix 2. The integration over intermediate variables yields a drastic simplification due to the appearance 
of quite a few (5-functions. This tedious but straightforward derivation is given in detail in Appendix 3. 



Starting with the general formula (A. 25) obtained in Appendix 3, we first evaluate the high energy behaviour of the 
integral over ha,b the dominant contribution of which, as well-known from the ordinary BFKL analysis, is obtained 
for Uafi = 0. The integrals are dominated by the usual BFKL saddle-points at ha^ ~ 5; ijj(^) = 4 log 2. Taking into 



3 



account in the saddle-point derivation the logarithmic corrections due to the zeroes in the prefactors at ha^b = 5, 
gets: 



Formula ( [lO| ) is then reported in the full amplitude (H). 
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(10) 



A few comments on our result ® are in order: 

i) Conformal invariance is explicitly realized, due to the decomposition on the combination of SL{2, C) Eigenvectors 



Eq{p')Eq{p). This property being shared by the tree-level amplitude, see (A^), is thus valid for the whole amplitude 
including one-loop level. Hence, saturation which is expected from the compensation between tree level and 
one-loop due to the different energy dependence is genuinely a conformal invariant phenomenon. We conjecture that 
this will be true to all orders, since only tensorial SL(2,C) properties play a role in the symmetry properties of the 
amplitude. Thus the saturation regime based on dipole loop contribution, if mathematically convergent, is expected 
to respect the conformal symmetry. 

ii) The key difference between the tree-level and one-loop amplitude properties lies in the difference in the analytic 
singularities of the integrands. While the tree level term in (^) is dominated by a moving saddle-point at high 
energy, the one loop amplitude ( [l0| ) has double poles at x{h)— 2x(l/2) =8 log 2 which depend only on the anomalous 
dimensions h and thus are independent of the energy. This very peculiar property means that the one-loop contribution 
may depend (and thus can give information) on BFKL properties which are not reached in the conventional analysis. 
In particular, all higher conformal spins may contribute, contrary to the tree level contribution. Also, contributions for 
values of h not limited to the interval ^] can give rise to poles. This may be interpreted as unlimited contributions 
of higher twist terms of the BFKL solution. 

iii) Apart the factor a'*, which is the natural strength scale of the one- loop amplitude, there is a dynamical factor 
connected to the square of the function 937? (/i, 5, ■^). This function is the QCD analogue of the triple Pomcron coupling. 
However, it is not the value gzvO^i \i \) which is here relevant, as was the case in the QCD dipole analysis of the 
triple pomeron coupling in |25|, since /i = i is not solution of the double pole equation. One could also note that the 



square of g^-p appears, contrary to saturation equations on an extended target, where the dependence is linear [hTl 



III. SATURATION CRITICAL LINE 



As discussed in section I, saturation is expected to take place when the one-loop amplitude, with coupling of order 
a^, becomes comparable in strength with the tree level amplitude of order . Formula (^) is completely explicit 
and thus ca n ser ve for analyzing saturation depending on transverse momentum or, equivalently using the Fourier 
transforms (A.4), on impact parameter. For sake of simplicity we will concentrate on the forward contribution at 
(5 = 0, that is (through unitarity) for the total cross-section integrated over impact parameter. In this limit, the 
SL{2,C) Eigenvectors simplify |23] to give 



E, 



P 



(11) 



Consequently, we may infer that the key "order parameter" for the discussion of saturation will be the scale invariant 
ratio of incident dipole sizes p/ p' , in conformity with the conformal invariance of the amplitude. This will in turn be 
determined by the ratio of the characteristic scales of the incoming systems which are encoded in their corresponding 
QCD dipole wave functions, see (||). 

The general solution of (^ is easily obtained by a summation over the double pole contributions verifying x(^i)^ 
2x(^), whose location are caracterized by an infinite number of values hi) = {ni,ji). One finds a formal double series 



(-1)^ 



[-"my 



(9l0g|ff3-p(/j., \, \) 

dh 



(12) 



h=hi 
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where, for sake of simplicity, we used the notation L = log \pl p'\ ■ 

For incident dipoles of comparable size, L ~ we may note that the one loop result then depends on an 

infinite sum over (71^,7^) which is due to the analytic properties of the function x(^)- It even raises the question 
of the convergence of the double series (^). On a physical ground, since the anomalous dimensions hi have no 
obvious limitation, this situation can be interpreted as due to a dependence of the one loop amplitude on all higher 
twist contributions, including those with non zero conformal spins n^. This is a quite strange example where higher 
twist calculations are required to all orders in a purely perturbative framework, in striking analogy with actual non 
perturbative situations. 

Leaving this regime for further work, we will here focus on the case where only the first pole with h* = (7*, n* = 0) 
contributes, i.e. where the equivalent of a leading twist contribution is sufficient. This approximation is realized in a 
region where L ^ 1 (or L <^ 1 by symmetry) . 

The dominant double pole in ( p^ such that x{h*) — '2x{\) — 8 log 2 is h* = {n* = 0,7* .183). Consequently, 
factorizing the tree level coupling, the resulting one loop amplitude can be written 



fSUL,Y) ^ ^ \pp'\ exp 



2Y + 2{h* 



Yi -31ogy + logi 



(13) 



where we define Y = uj{^)Y, and 



Y, = log ■ 



16^2 



53p(/i*, 2' 2) 



XC2) 



(14) 



which is a priori large, being in the perturbative domain a <C 1 and constant up to 1/L correction terms. 

We define the saturation critical line as the line in the {Y, L) where the tree-level contribution equals the one loop 
one. It is a simple extension to the conformal invariant configuration of the usual one 0]. Beyond this line, saturation 
is expected to take place. 

We now make use of the well-known saddle-point approximation of the tree level amplitude 



f'iU{L,Y)^^\pp'\ exp. 



^^^^\y + 2[h,-l/2)L~\og 



1 



x(l/2)' 



- 1/2 logy 



(15) 



where the tree level saddle point he is given by the well known implicit equation given in (^. Inserting the results 
( p^ ) and (jlj) in (||) with a redefinition of the constant 



= Fi + log 



1 



'27rx(i) 



we get the equation (O) announced in section I for the saturation critical line. 



(16) 



IV. DISCUSSION, OUTLOOK 

Let us describe general features of the critical line, when a is very small but otherwise not specified. Following 
expression (^, it is straightforward to realize that there is an absolute bound when L, F 3> 1 given by Y /2L > 
x{h)/x'{h*) ~ .09. Indeed below this bound, the compensation can no more take place between tree level and one 
loop contributions. It is interesting to note that the bound implies that a renormalization group evolution, with L 
moving to large values while keeping Y fixed, cannnot lead from a weak density to a saturation region, as expected 
from the general theoretical picture of the QCD dynamics of partons. 

One can delimit three typical regimes in the saturation transition, depending upon the behaviour of the tree level 
amplitude: 

I -C Lo— > Li <C : in this region, saturation is obtained in the region y > + 2£ — /i*) ^ Fo + .63L, 
where Yq is a large effective rapidity scale incorporating the logarithmic prefactors. 

II Lq <C Li— > ^2 < : It is well-known that the BFKL behaviour discussed in I acquires an extra "diffusion 
term" (of the gluon transverse momenta along the BFKL ladder contributions) in the intermediate region. Within 
the saddle-point approximation it leads to the following saturation curve: 
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1 u*\ X(^) 2L 



2 



r>vi + 2.(--rJ-^_, (17, 



where the last term is due to diffusion. Quite interestingly, saturation is favoured by diffusion. Indeed, formula (|1 
approximately interpolates between the behaviour I and II by a smooth decrease of the limiting slope of Y versus L 
III L2 ^ Y < L : near the bound 2L/Y < x{^)/x'{h*) ^ 10. In this region, one can use the relations 



Using this and again introducing the effective large rapidity scale Yq one obtains, within a good approximation at 
large (Y.L), the constraint 




y><^ -TTT + \^} ■ (18) 



It is easy to check that (|l|) satisfies the bound Y /2L > x (|) /x'{h*). 

We have not yet discussed the role of the logarithmic prefactors. Using our explicit expressions, it will be easy 
to incorporate them in a more complete study which we leave for the future. They may play a phenomenological 
role. Indeed, it is important to note that the lar ge fac tor | \ogY of (0) is due to the double zeroes appearing in the 



loop couplings to the intermediate dipoles, see ( |A.25| ). This has the important effect of a strong energy dependent 
damping on the effective triple Pomeron couplings, which could otherwise be large ]25[ |. Another possible important 
role would be, by mere extrapolation to higher dipole loops, to establish a useful hierarchy between these different 
orders. Indeed the coupling strength expected at order n, namely ~ (a^ exp a;(^)F)"' cannot provide a hierarchical 
structure in the saturation regime. We are thus conjecturing a series of critical saturation lines corresponding to 
the onset of compensation between a growing number of higher loop order contributions. The study of this regime, 
necessary to understand the full saturation behaviour in the conformal invariant set-up will be pursued in the near 
future. 

In fact, our work suggests a strategy to attack the full problem of saturation from dipole loops. It seems possible 
to analyze the n-loop dipole diagrams in the same way as we have shown here for the 2-loop case. It is clear that 
many diagrams can contribute, depending on the multiple rapidities where new dipoles can branch. However, we 
may conjecture that the integration over these intermediate rapidity intervals would dominantly lead to the QCD 
multi-Pomeron exchange and the corresponding 1 ^ n dipole vertices which were analyzed in Ref. | |28[ |. Then, it will 
be possible to explore the conformal invariant saturation regime it self and not only the transition to it. 

Another point merits some care, when the onium typical sizes are of the same order. A quite interesting regime 
seems to set in when the sum (|2|) over all the double poles in the one- loop integrals becomes relevant. While being 
still at small coupling constant, this regime shares with the non-perturbative one the features of a summation over 
higher twist contributions. It would be nice to solve this case in detail. 

Finally, the step between our calculations and actual experimental possibilities, which is not done in the present 
work, could be done either through factorization properties or by the consideration of "hard" subprocesses with very 
high energies like Mueller-Navelet jets at LHC or 7*— 7* collisions at a future linear accelerator. Relevant studies are 
certainly deserved in the future. 
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APPENDIX 1 



SL{2, C)-invariant formalism 

Using the 5i(2, C)-invariant formalism, the solution of the BFKL equation obtained at leading log level of pertur- 
bative QCD [|ol reads 



/q°\p,p'\Y) = ^ \pp'\ I dh EQip')E^{p) d{h) 



ih)Y 



(A.l) 



where the SU ( 3) co uplingct is fixed, small but unspecified, in the leading log approximation we consider. 

In equation (A.l), the symbolic notation J dh = X]^^-oo / '^'^ corresponds to the integration over the quantum 
numbers associated to the continuous unitary irreducible representations of SL(2,C), namely 



1— 1 + n ^ 
h = lu ^ : n = lu ^ : n £ Jv . £ Tc 



(A.2) 



E'q{p) and ujlh) are, respectively, the SL{2,C) Eigenfunctions and Eigenvalues of the BFKL kernel |23|. The Eigen- 
values read 



where * = (logE)'. The SL{2,C) Eigenvectors are definedQby 



P , P\ 



with 



(A.3) 



(A.4) 



p\ _ \p\ m 



d'r e-'Q^ EUp) - (-) 



h-h 



where h ^ 1 ~ h, r is the 2-d impact-parameter, and introducing useful notations 



d{h) 



.2 I (n-iy 



,2 I ("+!)' 



; b{h) 



\n\ 



a{h) 



27r2 



,2 I lil 



(A.5) 



(A.6) 



We will also make use of the formalism in the full coordinate space, the 4-diple amplitude T'^^^ being defined at 
tree level of QCD dipole contributions by 



r(°)(po,Pi;p[),p'i|r) ^^Jd'Q « fS\p,p'\y) 



(A.7) 



where the impact para mete r is r = h{Po + p'l^ Po — Pi)- 



Inserting (A.l) into (A.7), one finds 

TW(po,Pi;p^,p'i|r) = ^ Jdh G'^{po,pi;p'o,p[)(^^' + ^^ d(ft)e-('^)^. 



(A., 



where , using the 5i(2,C) Eigenvectors in coordinate space ( |A.5| ), the Green function G'' in coordinate space 23 
can be written as 



'^Note that an analytic expression of the Eigenvectors Eq{p) in the mixed representation has been provide d [Blj in terms of 
a combination of products of two Bessel functions. For s implicity, we did not include the impact factors pslpiJ^Phenomeno- 
logically, the leading contribution to the amplitude (A.l) comes from the n = component which corresponds to the BFKL 
Pomeron. 
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G\p^,p,-p'^,p\) ^ I d'po.E\po,a,Pi'a) E\p^^,p,^) = \pp'\ I d'Qe-^'Q eI{p')E^q{p) . (A.9) 

We will also define for convenience the "instantaneous" (i.e. Y=0) amplitude 

r(po,Pi;p;„p'i) = T(o)(po,Pi;p[„p'i|i" = 0) . (A.io) 

Due to completeness and frame independence of the formalism, many formulae can be simplified using these "instan- 
taneous" amplitudes between dipoles. 

The single multiplicity distribution of dipoles ni after an evolution characterized by the rapidity Y is simply related 
pl[ to the previous expressions (e.g. ([A.9|)), namely 



n,{p^,p,-p',,p[\Y)^^^ [ dh G"(po,Pi;Po,p'i)(^^'+"V4)e"('^'' ■ (A.U) 



Let us finally note a useful completeness formula valid at tree level: 

^(0) f ' ' IV- _L '\ - 1 f d^Paod^Pa, d^Pa'^d^Pa[ 

' {poPi;PoPi\Y = y+ya+y) ' 



X "1 (P0Pi;PaoPai|y) "1 (PoPi;PaoPaily') T^-^^Pao Pai , Pa'^Pa^lVa) , (A.12) 

where po,pi are the transverse coordinates of one initial dipole (resp. PoPi for the second one), PaoiPai are those 
of the interacting dipole emerging from the dipole po,Pi after evolution in rapidity y (resp. Pa'giPa[ and y' for the 
second one), which interaction is described by the amplitude T{pao, Pai'-, Pa' , Pa'^\ya) as in ( [A.S| ). 

Equation ( A. 12 ) has the physical interpretation that, at tree dipole level, the multiplication of dipoles due to 
rapidity evolution does not lead to saturation, but to the same tree level amplitude in a kind of bootstrap property. 
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APPENDIX 2 



The two-dipole density distribution n2 



The two-dipole density distribution n2 results from the solution of a sp ecific evolution equation which consists of 
an extension of the mixed evolution of the one formulated in ref. P, |3l| , p^ , namely: 



"2 {POPI, PaoPai, PboPbAv~V,v) 



aNc f dhdhadhf, 
2^ 



\PaPb\ 



duji- 



i(y-v) 



dhJ 



UJi — LU (h) J LO [ha) + Ll> (/if,) — UJ 



(A.13) 



where, by convention, diVi du)i/2i: and 



n2 (poPi;PaoPai,PhoPbi) = -777— TTTTTTTT I d'^pad^pii d^pj 7^|^;'^°^'" 

a{nja(na)a[nbj J 

) E^^{pb,(3,pb,p) E''{poj,pi^) 



with 



7^' 



h.h^M - f d^rad^rid^r2 



a, ^,7 



^02 7*12 



(A. 14) 



(A.15) 



where p = PQ-pi,Pa Pag-Pai, Pb Pbo-Pbi- 

Conformal invariance implies the well-known tensorial reduction |2£ 



rrjh^ha-hh r ^h—ha—hh r 



[Pia\ 



r_ iha—hij—h r_ -\hiy—har-h n i i \ 

LP/37J [Pja\ 93V (n,ha,hb) , 



(A.16) 



where g^-p {h, ha, hb) happens to be [ p5[ the QCD triple Pomeron coupling as obtained in the dipole formulation, 
namely: 



93V {h,ha,hb) 



d^rod^rid^r2 , 

\ ?2 [^01] 

Foi ro2 ri2\ 





ha r 


?*02 




To ^2. 





'*12 



(l-ri)(l-r2) 







ha r 




rQ2 












J0^2. 





^12 



(l-ri)(l-r2) 



(A.17) 



Considering the Fourier transforms (A.4) for the SL{2,C) Eigenvectors E^ in the mixed representation, one writes 
nY'°^'^''{poPi;PaoPai,PboPbi)^93v{h,ha,hb) ^^^^^j^^^^^b) ^pp^p^^ J (fq^cfq^d'^QS^'^\Q+qa+qb) 

X E'^{p)El'^{pa)E^^^{pb) e-*('?»'-»+9''''''+'3'-) J d^v d^w e-<(«°^^)^+Of ) 
X [(y;_„)/2]-i+''-''''-+^'' [{w + i;)/2]-i-^+''°-'^^} X {a.h.} , (A.18) 

where p^p — v, pa^ + pp-y — w, 2ra = Pao + Pai, 2rb — pbo + pi and r is the overall impact parameter. The notation 
{a.h.} indicates the anti-holomorphic part of the bracketed term in the integrand for which the integration variables 
are complex conjugates and the expo nents hi are replaced by hi. Equivalently, the distribution n2 for the lower vertex 
is given by the same equation ( A.18 ) by using prime indices. 
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APPENDIX 3 



The general one-loop amplitude 



Implementing equations( A.8) and (A. 18) in the master equation and first integrating over impact parameters, 
one gets: 

(27r)8 5{qa - q\)d{qt - q\)S{q', - q\)5{q', - q" t) . 
Then, integrating over the intermediate dipole sizes, and using the orthogonahty [p3|| relation 



where e"^^''^ = Ej^ (p) / E^'^^ (p) is just a p, g- independent phase factor given in Rcf. |2^, one finds 

(2^)« 24 <S(2) (ha , h\)6^-^^ ih,h\)d^'\hll-h\)S^^\hll~ /l" ,) , 

where 5^'^\h^h') = Snn'^iv — v'). The integration over qa — qt finall y gives {2tt)^ 6{v — v'). 

Plugging in formulae (^,^) the result obtained in formula (A. 18) and integrating over (5-fuctions, one gets: 



(A.19) 



{p,p'\Y = y + y') = - 
X / dhdh 



dhadhb 



d{ha) d{hb) 



2n^J J HKW HhbW 

, ff3P [h, hg, hb) gl-p{h' , hg, hb) , „h / ^ -^h' , 



Y-y' 



h{l-h) 

uy+ui{y~y) 



h{h') 



pp'\E^{p) Eq{p') TiQ^hX) 

^yy'+uj[{y'-y') 



X I dy dy' / dujdujiduj'du!^ -, — 7- — ^ 7- — ^ t-, tttt -, — v, — ^ v, — ^ — ; , 



(A.20) 



where, after the change of variables w — v{\ — 2i), w' — v(\ — 2t'), 
nQih,h') ^ f d'^vd'^tdh' e'Q"'-*-''^ i^yh-h'-i^-Wi-iii,-h,,^-^^_^yWi-h,-i^,,^,h^-h,-h'^^ ^ ^ (^ 21) 



After integration over v one obtains |30 

h—h' / n \ h — h 



■HQih, h') = e^^("-"') ^{h - h') I dh dh' (t t') 



l\h' -h 



X {a.h.} 



(A.22) 



where, by definition, 7(2) — r(z)/r(l — 2) . 

The remaining integral can be proven]^ to reduce to 



Hq(/i, h') 



2 \ h—h' / 2 ^ h—h' 



(2/1-1) (2/1-1) 



[S{h-hf) + i-)"^^+''' -fih + hg-hb) ^{l-h'-hg + hb) Sil-h-h')] . 

(A.23) 



The integral over h' can now easily be performed. The remaining poles a,t h' — h and h' = 1 — h give twice the same 
contribution due to the over completeness relation of the Eq generators One finally obtains: 



^One shows using general methods that TlQ{h, h') is a combination of 5 functions whose coefficients can then be computed 
directly on formula (A.22). 
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27r2 



dhdhadhb 



g3v{h,ha,hb) 



dy 



dy' / dujduiiduj' dui'i 



b{ha)h{hb)h{h) 



d{ha)d{hb) 



pp'\ E^ip) E^ip') 



^u)'y' +uj[{y' -y') 



{uj{ha)+Uj{hb)-Uj){uJi-Uj{h)) {uj{ha)+LLl{hb)-U!'){uj[-Uj{h)) 



(A.24) 



This one-loop amplitude preserves the global conformal invariance of the tree-level BFKL 4-gluon amplitude, since 
the only scale-dependence on Q is present in the conformal Eigenvectors \pp'\ Eq{p) Eq{p'). 

It was noted in Ref. |2|] that the integration over rapidity variables yields two different contributions depending 
on the sign of the quantity uj{ha)+uj{hb) —u}{h). Indeed for u}{ha)+u!{hb) <uj{h), the relevant poles are situated at 
uj = LJi ^ uj' ~ Lo'i — uj{h), leading to a contribution associated with the single Pomeron dependence e"^'*^ ^ . In the 



opposite case, namely Lu{ha)+(-o{hb) >uj{h), the relevant poles are situated at lj = uji - 
resulting amplitude 



-Lo' = uj'i = Lu{ha)+i-u{hb). The 



f^" {P,P'\Y) 



1 

dhadhb 



PP 



dh 



-ir e'q{p')e^{p) 



-1 2 



ha{l-ha)hb{l-hb) 



Ju>{h^)+uj{hi))Y 



g3-p{h,ha,hb) 



X(h)-x{ha)-x{hb) 



(A.25) 



is the one we are interested in for the saturation problem since it possesses the QCD "double Pomeron" energy 
behaviour ^i^iha,)M^{K))Y ^ ^i^ih) y ^ This is precisely the energy behaviour which will compensate for the higher order 
in the coupling constant. Notice that the expression depends only on the sum Y = y + y\ as it should do from 
longitudinal boost invariance. 
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FIGURES 




Figure 1 

One-loop dipole amplttude in transverse coordinate space 

The transverse coordinates pi of the incident and of the two 

interacting dipoles are simply denoted by their respective indices i. 

The coordinates 2, 2' refer to the integration points, see Appendix 

2. 




Figure 2 

One-loop dipole amplitude in rapidity space 

The initial dipoles poPi , PoPi evolve in the rapidity range y — 
y',y' — y' with one-Pomeron type of evolution and then with two- 
Pomcron type of evolution over the range y, y' . Note that y-\-y' = 
Y, is the total rapidity range. After integration over y,y', see 
Appendix 3, the resulting amplitude depends on the total range 
Y. 
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Saturation 



Y 




I 



II 



III 







L 



L 1 



L 2 



Figure 3 

Critical Saturation line 

The critical saturation line is represented in the plane Y , L. The saturation region is under 
the curve. The region beyond the vahdity of leading twist approximation is below a value 
1 << Lq. The regions I, II, III correspond to different regimes of compensation between 
tree and one- loop dipole contributions (see text). 
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